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1.0 SUMMARY 


An integrated autopilot/autothrottle system was designed using a total energy control 
design philosophy. This system was designed around a fixed inner loop configuration that 
controls flightpath angle (FPA) and velocity using elevator and thrust. Each vertical control 
mode: 

1) Glideslope control 

2) Vertical path control 

3) Altitude control 

4) FPA and go-around control 

5) Speed profile control 

6) Mach control 

7) CAS control 

8) Angle-of-attack control 

9) Groundspeed control 

10) Flare mode 

11) Velocity vector control column steering (vel-CCS) 

generates an FPA and velocity command that are compatible with the inner loop. 

The system design uses a “total energy control concept” which ensures that the system can 
differentiate between maneuvers requiring a change in thrust to accomplish a net energy 
change (e.g., an increase in velocity or height) and those maneuvers that only require elevator 
control to redistribute energy (e.g., speed/altitude exchange). 

Previous work under Boeing and NASA funding developed the linear design of the inner loop 
configuration and basic outer loops. From this linear model, a nonlinear simulation model was 
developed and the system was evaluated on a flight simulator. 

This report covers the further development of the system, in particular: 

1) Engine EPR controller redesign 

2) Flare mode and go-around 

3) Vel-CCS 

A summary of simulation results have been included to illustrate the performance of the 
various modes. 

The overall system design has been achieved with reduced complexity as compared with con- 
ventional autopilot/autothrottle systems, and yet is more versatile and achieves better 
performance. 


All modes of operation (except go-around) satisfy the basic performance criteria: 

1) Less than 1-kn speed error for any flightpath maneuver 

2) Less than 20-ft height error for any speed maneuver 

3) Not greater than O.lg normal acceleration (except go-around) 

The system was designed for maximum safety. In particular, the angle-of-attack mode was 
designed to prevent stall in all conditions including full extension or retraction of flaps. 
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2.0 INTRODUCTION 


The objective of the total energy control system (TECS) design was to develop all longitu- 
dinal control modes for the autopilot/autothrottle and flight management system using a com- 
mon, generalized flightpath and speed control algorithm, thereby eliminating unnecessary 
control law replication and providing a flight management system (FMS) simplification. 

This algorithm provides full-time thrust and elevator control coordination for all modes and 
flight conditions, thus eliminating speed deviation due to flightpath maneuvers and flightpath 
deviations due to speed maneuvers. 

The TECS was designed to provide a full complement of longitudinal FMS modes, including a 
properly prioritized and integrated stall, and engine overboost protection. These provisions 
have been designed to be foolproof. A maximum thrust stall, as happened on a Mexicana 
DC-10 during climbout, would not be possible. In addition, speed control has been prioritized 
during descent to prevent overspeeding when the throttle goes to idle. Present autothrottles 
cannot provide these operational safeguards. 

A significant aim of this work was to develop a total longitudinal system that takes care of all 
control requirements simultaneously, thereby providing a significant reduction in software 
and hardware coupled with an improvement in overall performance. 

The total energy control concept requires that thrust be used to control the total energy re- 
quirements of the aircraft while the elevator is used to distribute the energy between the 
flightpath and speed objectives. 

The system design approach is to: 

1) Develop an energy rate error (E f ) which is the sum of longitudinal acceleration error 
(V e ) and FPA error (y e ) (i.e., E c = V e + gy e where g = gravitational constant). 

2) Develop an energy rate distribution error (E 5 ) which is the difference between V e 
and 7 e (i.e., E 6 = V e - gy e ). 

The energy rate error is used to develop thrust commands, whereas the energy rate distribu- 
tion error is used to develop elevator commands. The commanded values y c and V c are devel- 
oped by normalization of the particular flightpath and speed control signals for each 
individual mode. 

Control of thrust based on total energy error also means that maneuvers involving mainly 
energy exchange are executed by the elevator and unnecessary throttle activity are 
eliminated. 
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This report documents the development of TECS and covers: 

1) A redesign of the engine controller to reduce steady -state gain and bandwidth varia- 
tion over a wide range of aerodynamic conditions and to prevent engagement 
transients. 

2) Design of a flare law and go-around mode. 

3) Design of vel-CCS. 

4) A performance evaluation to illustrate the current system capabilities and highlight 
the unique features of TECS (including the improvements to existing modes: (a) 
altitude select mode, (b) velocity select mode, and (c) AOA mode, and demonstrate 
the prioritizing of speed in double maneuvers that cause limiting of the throttle. 
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3.0 SYMBOLS AND ABBREVIATIONS 


AOA 

ACSL 

ALFREF 

ALTCMD 

ALW 

AZ, AZCG, A z 
CAS 

CAS SEL 

CCS 

CG, eg 
CMD 
DCOL 
DELCAS 
DEPRO 
DFM 

E e 

EAS 

ELEV 

EPR 

EPRo EPRcmd 
EPRic 
EPRidl 
EPRqut 


Angle of attack 

Advanced continuous simulation language 
Angle-of-attack reference value 
Commanded altitude 
Wing angle of attack 

Longitudinal acceleration (through center of gravity) 

Calibrated airspeed 

Selected CAS 

Control column steering 

Center of gravity 

Commanded signal 

Column deflection 

Error in CAS 

Input to engine control loop 
Flap angle 

Energy rate distribution error 
Energy rate error 
Equivalent airspeed 
Elevator angle 
Engine pressure ratio 
Commanded EPR 
Initial EPR 
Idle EPR 
Measured EPR 
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3.0 SYMBOLS AND ABBREVIATIONS (Continued) 


FL 

FLEPRE 

FMS 

F N C 

FPA 

FPA L1,L2 

G 

G l 

g 

GAINT 

GAINV 

GAMDED 

GAMEPS 

GAMMA 

GAMMAD 

GEAR 

Gis) 

GSE, GSERR 
H, h 
h 
h 

^B’ ^baro 
h c> h CMD 

h dot 


Flap angle 

Go-around engagement indication (1 = engage) 
Flight management system 
Net thrust command 
Flightpath angle 

FPA for leg 1 of maneuver, leg 2 of maneuver 
Engine gain (throttle command to EPR) 

Engine gain (XSA to EPR) 

Gravitational constant 
Gain within engine loop, identical to Kp 
7CMD f° r unit stick deflection 
Commanded flightpath angle (deg) 

Flightpath angle error 
Flightpath angle 
Flightpath angle (in deg) 

Landing gear position (1 = down, 0 = up) 

Forward path transfer function 

Glideslope error 

Height 

Height rate 

Vertical acceleration 

Barometric height 

Commanded height 

Height rate 

♦ 
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3.0 SYMBOLS AND ABBREVIATIONS (Continued) 


Hdtcmd 

^dtcrnl 

■^dtrl 

H ERR> h error 
hi 

h 0 

h 0 Ll 

h R 

href 

h S ei 

h TD 

h £ 

ICTIC 

KAEPR 

KALFBS 

K P 

MACSEL 

MCP 

RADALT 

RAND1 

RL 

SD 

SX 


Commanded height rate 
Limited commanded height rate 
Rate limited height rate 
Height error 
Inertial height 
Initial height 

Initial height for legj of flight 

Radar altimeter measurement, height above runway 

Reference height 

Selected height 

Height rate at touchdown 

Height error 

Engine loop integrator initial condition 
Controller gain equivalent to G 
Bias for computation of K a 

Forward path engine controller gain identical to GAINT 

Selected Mach number 

Mode control panel 

Radio altimeter measurements 

Range rate 

Rate limiter 

Standard deviation 

Touchdown along runway 
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3.0 SYMBOLS AND ABBREVIATIONS (Continued) 

T Time 

t Time, time from start of flare 

TAS True airspeed 

TCAB1 Hardware throttle position 

TECS Tbtal energy control system 

TH,THPOS,THR Throttle position 

THBIAS Throttle bias 

THCM Throttle command 

THETAD Attitude angle 

THMAX Maximum throttle angle 

THMIN Minimum throttle angle 

THTRM Throttle trim position 

TOTT Ibtal air temperature 

u g Longitudinal gust velocity 

Vcas> Vcasel Calibrated airspeed, selected CAS 

V CMD> Vcmd Commanded velocity, velocity rate 

VcmdFL Commanded velocity rate during flare 

V d Derived velocity rate 

Vel-CCS Velocity vector control column steering 

V err , V error Velocity error 

Vq> Vqq Groundspeed, initial 

Vq! Smoothed groundspeed 

V HF High-frequency component of velocity rate 

V LF Low-frequency component of velocity rate 
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3.0 SYMBOLS AND ABBREVIATIONS (Continued) 


V P 

V ref 

V T 

V £ 

w 

w 

g 

X 

XSA, XSHARG 

a 

a ref 

a vane> a w 
' Yc ’ Tcmd 


Tl 

7SYNT 


7e 


Ah G/s 

AM C 

AT 

AV C 

A 7SEL 

8 

S COL 

6e 


True airspeed 
Reference velocity 
True airspeed 
Error in velocity rate 
Weight 

Vertical gust velocity 

Distance along runway from glideslope intersection point 

Cross-shaft angle 

Angle of attack 

Reference angle of attack 

Wing angle of attack 

Commanded FPA 

Inertial FPA 

Synthesized FPA 

FPA Error 

Height error off glideslope 

Change in commanded Mach number 

Change in throttle position 

Change in commanded velocity 

Change in selected FPA 

Pressure ratio 

Column deflection 

Change in elevator angle 

Change in commanded elevator angle 
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3.0 SYMBOLS AND ABBREVIATIONS (Continued) 


r 

e 


%’ CT5 th 


r 

T E 


r F 


T GA 


Ty 


Energy 

Damping ratio 
Attitude angle 

Standard deviation of height rate, velocity 

Standard deviation of elevator angle, throttle position 

Time constant of engine loop, or steady state CCS lag 

Time constant of engine 

Time constant of flare 

Time constant of go-around filter 

Initial system time constant 

Time constant of V complementary filter 

Natural frequency 
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4.0 DESIGN AND PERFORMANCE REQUIREMENTS 


1) Integrate all longitudinal and vertical control modes to use a single generalized flight- 
path and speed control algorithm that provides simultaneous thrust and elevator 
commands. 

2) Provide full-time elevator and throttle control coordination to limit speed coupling errors 
due to flightpath maneuvers to 1 kn and limit flightpath coupling errors due to speed 
maneuvers to 20 ft. 

3) Eliminate throttle response to flightpath perturbations. 

4) Minimize throttle response in cases of desired altitude/speed trades. 

5) Provide capability to limit normal acceleration during vertical maneuvers to preselected 
value (i.e., 0.1 g). 

6) For step command inputs the response should be smooth and overshoot free. The damp- 
ing ratio of dominant poles should be greater than 0.7. There should be a consistent 
response over the whole aerodynamic envelope. 

7) Control mode switching should be transient free. 

8) The design should incorporate foolproof safeguards to prevent stall and overboost of the 
engine for all modes and flight conditions. 

9) For final approach, the system should be optimized for flightpath and speed tracking, and 
commensurate with acceptable throttle/elevator activity in turbulence. 

10) For cruise, the system should be optimized for minimal throttle/elevator activity and 
commensurate with acceptable speed and path tracking in turbulence and windshear. 

11) The flare mode should achieve an overall longitudinal touchdown dispersion <j x < 125 ft 
and a sinkrate dispersion 0 ^ < 1 ft/s, while exhibiting no appreciable sensitivity to ap- 
proach terrain and runway slope. 

12) For the velocity vector steering mode the requirements are that: 

a) The response lag between the pilot input and the displayed control variable (through 
which the pilot closes the short-term control loop) must be small enough to allow 
smooth and positive control. 

b) With the pilot out of the loop, the automatic control must track the control command 
close enough to avoid pilot intervention for the purposes of short-term airplane 
stabilization. 

c) During manual maneuvering, the automatic system must not in any way interfere 
with the pilot’s control task, but rather should provide the necessary command and 
stability augmentation to yield good handling qualities reflected in control column 
sensitivity, responsiveness, and damping. 
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5.0 ENGINE CONTROL 


5.1 ENGINE CONTROL OBJECTIVES 

Previous autothrottle systems have employed a simple throttle lever position servo control. 
The problems with this type of control were: 1) significant backlash between the throttles and 
the engine that resulted in limit cycling, and 2) widely varying engine static gain and re- 
sponse over the flight envelope which resulted in poor system performance. 

The TECS computes the total net thrust required to execute both flightpath and speed com- 
mands, basically for a nominal weight aircraft. The new thrust command is then scaled for the 
actual aircraft weight, which is presumed to be available from the flight management com- 
puter. To avoid the problems described above, the engine should be controlled to deliver di- 
rectly the total required net thrust. However, the available engine control parameter is engine 
pressure ratio (EPR). Therefore, the new incremental thrust command must first be converted 
to an EPR command (EPR C ). This is done by conversion of net thrust command Fn c to normal- 
ized net thrust command FNq/6 using the dividing factor 8 - p/p 0 . 

Next, the normalized net thrust command FNq/ 6 is converted to EPR C using the gain KEPR. 
This gain is derived from Figure 1, which shows the relationship between EPR and F N /6 for 
the JT8D engine adjusted for dual engines. This FNq/8 EPR relationship is Mach dependent; 
nevertheless, the Mach effect is small for this type of engine and the relationship may be 
approximated by an averaged constant. 

5.2 EARLIER ENGINE CONTROL USING FIXED-GAIN EPR LOOP 

In the earlier TECS engine control design, (ref. 1), control of net thrust was achieved by 
conversion of the thrust command into a corresponding EPR command with the use of a fixed- 
gain EPR feedback loop around the basic throttle position loop (fig. 2). Assuming zero throttle 
position corresponds to the desired engine idle, the bias EPRIDL must be subtracted from the 
EPR feedback to obtain proper trim. This design approach reduced the effect of throttle to 
engine backlash and variation in engine dynamics. It was satisfactory for the low-speed ap- 
proach flight condition. However, tests on the flight simulator revealed other flight conditions 
where a mismatch existed between the computed throttle command and the throttle position 
(THTRM) calculated in trim mode. This caused transients at the start of simulation runs. 
Table 1 shows the calculated throttle command (THCM) and throttle position (THTRM) with 
Case 3 showing a 13-deg mismatch. In addition, engine characteristics were found to vary 
with flight conditions causing variation in loop bandwidth and steady-state gain between EPR 
command and EPR. Since this variation in steady-state gain was excessive, it was necessary 
to redesign the engine loop gains to reduce this variation. 


PRECEDING PAGE BLANK NOT FILMED 
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Table 1. Initial Throttle Command and Trim Value 


CASE 

NO. 

AERODYNAMICS 

CONDITIONS 

THROTTLE 

COMMAND 

THROTTLE 
TRIM VALUE 

EPR 

(initial 

condition) 

EPR 

(at 

idle) 


Altitude 

EAS 

THCM 

THTRM 

EPRIC 

EPRIDL 


(ft) 

(kn) 

(deg) 

(deg) 



1 

5,000 

150 

20.18 

21.02 

1.386 

1.025 

2 

10,000 

150 

25.1 

22.5 

1.484 

1.036 

3 

20,000 

150 

40.1 

27.7 

1.775 

1.004 

4 

20,000 

250 

20.5 

20.2 

1.436 

0.925 

5 

30,000 

250 

32.7 

24.4 

1.68 

0.97 


5.3 ANALYSIS OF SYSTEM 

Tb gain a better understanding of the problem, the system was analyzed using linear tech- 
niques. However, since the JT8D-9 turbofan engine has a highly nonlinear response caused by 
variable rate limiting during spool up and spool down, it should be noted that any meaningful 
conclusions drawn from the linear analysis must bear these factors in mind. 

5.3.1 First Order System 


An earlier analysis of TECS using a linear system assumed that the dominant response of the 
engine can be modeled using a simple lag with a time constant (t e ) of 1 sec. This approxima- 
tion was substantiated by studies which obtained favorable comparisons between linear simu- 
lation results of the full system and simulations using the flight simulator which employs the 
detailed JT8D-9 model, (ref. 1). 


It appears logical to start a linear analysis of the EPR command loop by assuming that the 
engine can be modeled using a lag with t e = 1 sec. In this case, the effect of the lag of the 
throttle servo will be ignored. 

For the system shown in Figure 2, the closed loop transfer function (CLTF) between EPR 
command and EPR is given by: 

AEPR KEPRP.G 

AEPR cmd = (1 + KEPRP.G.KEFB) + TeS 

where G is the engine gain, throttle position (deg) to EPR. By inspection, it can be seen that by 
defining 


K„ 

KEPRP = 

G 


( 2 ) 


Equation (1) can be simplified to: 

AEPR Kp 

AEPR cmd (1 + K p .KEFB) + t e $ 
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For this system, it is required that the steady-state gain (SSG) AEPR/AEPR CMD is 1. 
Therefore, 

SSG = = K P 

(s^O) 1 + Kp.KEFB 

Dividing top and bottom of the RHS of Equation (3) by 

(1 + K^KEFB) 

yields an expression for the time constant (r) of the loop 

t e 

T “ 1 + K p .KEFB 

Rearranging Equation (5): 


(4) 


(5) 


— = (1 + K KEFB) w 

T 

and substituting in Equation (4): 

K p = - (7) 

T 

from Equation (6) and rearranging gives an expression for the feedback gain: 

1 

KEFB = (1 - — ) 

or using Equation (7): 

KEFB = (1 - — ) 

t e 

Knowing the time constant of the engine lag, the forward and feedback gains of the EPR loop 
are dictated by the required overall time constant of the loop (i.e., equations (7) and (9). 


(8) 


(9) 


For an engine time constant (r^) of 1 sec: 

1 



Hence from Equations (2) and (9): 

KEPRP = 

Gt 

KEFB = (1 - r) 


(10) 


(ID 

( 12 ) 
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5.3.2 Second-Order System 


A study of the simulation JT8D-9 model used with the flight simulator showed that although 
rate limiting effects dominate, the linear response to very small perturbations is modeled by a 
lag of t e = 0.2 sec. This time constant is the same order of magnitude as the throttle servo, 
and therefore the effect of the servo should be considered in providing a comparison with the 
first order analysis. 


In addition, it is convenient to redefine the gain of the engine as: 

_ AEPR 
Gl = AXSA 

where XSA = cross-shaft angle. 


(13) 


The relationship between throttle (T) and XSA is approximated (in the linear region) by 

AXSA = 0.9 AT (14) 


Therefore, 


AEPR AXSA AEPR 
AT ” AT AXSA 


(15) 


This provides a more precise analysis and is convenient in that all figures given later in this 
section show engine gain (Gj) as AEPR/AXSA. 

The transfer function for this system is given by 


KEPRP • 0.9 • G x 
AEPR (1 + .15s)(1+t e s) 

AEPR cmd “ 1 + 0.9 • (KEPRP • G t • KEFB) 

(1+0.155)(1+t e 5) 


Redefining KEPRP = Kp/G 1( for the second order system analysis and multiplying 
by (1 +0.15s)(l + t e s). 


Then 


AEPR 0.9 K p 

AEPR cmd = (l+0.15s)(l+ t e s) + 0.9 K p KEFB 


(17) 


0.9 x K p /0.15r E 

^2 t e + 0.15s 1 + O^K^Kjjpg 

0.1 5t e 0.1 5t e 
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We require the SSG to be 1, therefore 


1 - 

KeFB 


Q-9K p 

1 + 0.9KJCEFB 


= (1 - 


0.9K„ 


(18) 

(19) 


The elimination of engine gain from Equation (17) and the use of Equations (18) and (19) 
shows that SSG can be maintained at unity over the whole flight regime provided the feed- 
back gain is selected as a function of forward path gain. 


The general equation for a second-order system is: 


s 2 + 2 


( 20 ) 


where 


oo n = natural frequency (bandwidth) 
f = damping ratio 

Comparison of this equation with Equation (17) yields 

t e + 0.15 


2 fan = 


0.1 5rr 


1 + 0.9K p KEFB 

0.1 5tp 


( 21 ) 

(22) 


Using Equations (21) and (22), it is possible to calculate the required forward and feedback 
gains K p , K EFB for selected values of damping and engine time constant t e to satisfy Equation 
(18). These results are shown in Table 2. 


Table 2. Results of Second-Order Analysis 


CASE 


T E 


K p 

K EFB 

LOOP GAIN 
0.9 K p K efb 

1 

1.0 

0.2 

5.8 

1.135 

0.02 

0.020 

2 

0.7 

0.2 

8.3 

2.316 

0.52 

1.08 

3 

1.0 

1.0 

3.8 

2.44 

0.55 

1.20 

4 

0.7 

1.0 

5.5 

5.0 

0.78 

3.50 
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Conversely, Table 2 shows that by selecting K p = 2.37 and K EFB = 0.52 (loop gain 1.23), the 
engine time constant can vary from 0.2 to 1.0 while damping will only change from 1 to 0.7. A 
5-to-l variation in engine time constant is reduced to a 2-to-l variation in overall loop 
bandwidth. 

5 . 3.3 Compensation for Variation of Engine Gain G 

A simplified diagram of the JT8D-9 engine is shown in Figure 3. The gain (G) across the 
engine (EPR/AT) is defined as: 


G = 0.9 X 


A(EPR) 

A(XSA) 


(23) 


The relationship between EPR and XSA is given in Figure 4. It can be seen that gain G is a 
function of total air temperature (TOTT), Mach number, and bleed valve position. Figure 5 
shows that the major factor in gain variation is TOTT. This dependency of G on TOTT has 
been plotted as gain factor KAEPR in Figure 6 for closed bleed valve and median Mach 
number as determined from Figure 4. 


Thus, to maintain unity steady-state gain of the transfer function 

AEPR 

aepr cmd 

and minimize loop bandwidth change with variation of the engine response time constant, the 
gains K p ( = KEPRP.G) and K EFB must be maintained constant. This requires that the forward 
path gain KEPRP be varied inversely with the variation of KAEPR. Thus, 


KEPRP = 


2.37 

KAEPR 


Further, it can be seen in Figure 4 that for small throttle angles 0 to 7 deg (or XSA between 38 
and 45 deg) the engine gain 


AEPR 

AXSA 

is very nonlinear. Therefore, to achieve satisfactory closed-loop performance, it was decided to 
use a lower throttle position limit of 7 deg and thus determine EPRIDL by taking the steady- 
state EPR at 45-deg cross-shaft angle from Figure 4. EPRIDL is programmed as a function of 
Mach number only. At the same time, a 7-deg bias is added to the throttle command to provide 
proper trim of the control loop. The resulting revised engine control loop is shown in Figure 7. 


Comparison of Figures 2 and 7 shows that DEPRO corresponds to <5EPR CMD in Figure 2. The 
output of the thrust integrator (“A”) corresponds to the nominal net incremental thrust com- 
mand computed by the TECS control algorithm. This signal is scaled by the airplane weight 
factor (W/W 0 ) and the atmospheric pressure ratio (1/5 = P(/p) to generate the normalized net 
thrust command (Fn c / 5), which is converted into the equivalent 5EPR CMD using the gain 


18 



KEPR. This gain KEPR is determined from Figure 1, using half the averaged slope for vary- 
ing Mach number. Half the averaged slope is used to account for two engines producing thrust. 

5.3.4 Thrust Control Command Initialization 


It is necessary to initialize the thrust integrator (“A”) to develop a AEPRqmd (DEPRO) corres- 
ponding to existing engine trim at engagement of the system. On the flight simulator, trim 
thrust and throttle position (THTRM) are calculated by the trim routine. This throttle position 
may be used to calculate the initial value ICTIC of the thrust integrator (“A”). 


From Figure 7 throttle command: 


THCM = THBIAS + KEPRP (DEPRO - KEFB (EPROUT - EPRIDL)) (24) 

During trim then: 

THCM = THTRM (25) 

Rearranging Equation (24): 

THTRM — THRTAS 

DEPRO = KEPRP + KEFB (EPR0UT " EPRIDL ) (26) 


From Figure 7 (during trim): 

DEPRO = ICTIC. KD.KEPR 

Therefore, 


ICTIC = 


DEPRO 

KD.KEPR 


(27) 

(28) 


Substituting in 

ICTIC = 
ICTIC = 


for DEPRO: 

(THTRM - THBIAS + KEFB (EPROUT - EPRIDL) KEPRP 

THTRM - THBIAS + KEFB (EPROUT - EPRIDL) KEPRP 

KEPRP.KD.KEPR 


(29) 


This provides a convenient method of calculating the integrator IC (ICTIC) utilizing the throt- 
tle trim value. 


5.3.5 Thrust Limiting 

To prevent engine overboost, the computed AEPR command is combined with the EPRj DL 
signal to form the total EPR command signal. This signal is used to form EPR error signals 
relative to the maximum and minimum allowable EPR values. If the EPR command exceeds 
the limit, the error signal is fed back through a high gain into the thrust integrator, thereby 
preventing the EPR command from exceeding the EPR limit significantly (fig 7). 

Likewise, the throttle position command is checked against upper and lower limits, and in the 
case of exceedance of the limits, the error is fed back through a high gain. 
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5.4 RESULTS OBTAINED FROM THE HARRIS SIMULATOR 


The linear analysis results shown in Table 2 (Case 2) were used as a starting point for the 
evaluation of the system on the flight simulator. This simulation included throttle servo and 
cable system nonlinearities. These nonlinearities caused a deterioration of the damping of the 
EPR control loop. To improve the damping, the loop gain 0.9K p .KEFB was reduced to 0.5, 
yielding for the forward path gain (0.9 K p = 1.5) and the feedback path gain (KEFB = 0.333). 
As a result 


KEPRP = 


1.5 

KAEPR 


where 


KAEPR = flTOTT) varied as shown in Figure 6. 
Table 3 details the results obtained from the flight simulator. 


Table 3. Results Obtained From Harris Simulator for 0.2 Change in EPR Command 


INITIAL VALUES 

STEADY-STATE VALUE FOLLOWING STEP IN EPR CMD 

CASE 

ALTITUDE 

(ft) 

CAS 

(kts) 

THCM 

(deg) 

THTRM 

(deg) 

EPRIC 

EPRIDL 

n 

DEPR0 

THCM 

(deg) 

XSHARG 

(deg) 

EPR0UT 

s.s.g. 

EPR 

EPRCMD 

d(EPR) 

imky 


1 

5000 

150 

21.02 

21.02 

1.39 

1.03 

57.8 

.365 

28.6 

65.2 

1.59 

1.0 

.0277 

.0249 

2 

10000 

150 

22,49 

22.49 

1.48 

1.04 

59.2 

.444 

29.5 

66.1 

1.69 

1.0 

.0296 

.0266 

3 

20000 

150 

27.75 

27.75 

1.78 

1.00 

64.3 

.711 

33.9 

70.1 

1.97 

0.95 

.8332 

.0299 

4 

20000 

250 

20.47 

20.47 

1.44 

0.93 

57.3 








5 

30000 

250 

24.37 

24.37 

1.68 

0.97 

61.0 

.694 

31.1 

67.7 

1.90 

1.1 

.0334 

.03 


Refer to Figure 7 for definition of variable names. 


It can be seen that the throttle command (THCM) exactly matches the trimmed throttle posi- 
tion (THTRM). For step responses over the range of aerodynamic conditions shown, then the 
steady-state gain variation was less than ± 10%. Figure 8 shows step response (AEPR cmd = 
0.2) for both the software and hardware simulations of the throttle. An excellent response was 
obtained using the software simulation, whereas the hardware caused minor oscillation. At- 
tempts have previously been made to provide an exact match between the software and hard- 
ware throttle by matching the rate limits and linear response. 
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Figure 9 shows the result of a ramp input of 0.04 EPR/s. Again a minor oscillation was evident 
with the hardware, though of reduced magnitude. As the input to the EPR control loop is via 
an integrator, then for most types of command changes to the system, the input to the EPR 
loop will take the form of a ramp. In practice, the small oscillation has not caused any notice- 
able problems. However, should this happen, then reduction of loop gain will improve the 
response. 

Results at 20,000 ft, 150 kn (Case 3) yielded similar performance (fig. 10). Figure 11 shows the 
performance at 1,500 ft, 120 kn, flaps 40-deg, and gear down for a series of step commands 
(DEPRO) designed to compare normal spool up, spool down, and spool up after sustained 
operation with the throttles at idle. Although spool up is considerably slower, EPR shows no 
overshoot. 


21 (Page 22 is blank) 


6.0 FLARE LAW 


6.1 INTRODUCTION 

Conventional flare laws in use today control sink rate as a function of altitude above the 
runway (ref. 2). A simple form of this law consists of summing radar altitude, sink rate, an 
altitude bias term, a ramp (to provide average trim requirements), and inner loop damping to 
generate an elevator command (<5e c ) during flare. For example, 

8e c = K h h + K h (h + h B ) + Damping + Ramp (30) 


where 


be c = elevator command 

h = vehicle altitude 

h B = bias signal 

— h = sink rate 

K h , K h and h B = designed gains 

Ignoring the trim and damping requirements, then the elevator is driven to satisfy: 

K „h + K h (h + h B ) = 0 (31) 

At decreasing altitudes, the control law calls for a proportional reduction in sink rate. The 
bias term h B provides the desired sink rate at touchdown. 

The solution of Equation (31) is: 

h{t) = ih 0 + h^'r- h B (32) 

where 



h Q = altitude at flare initiation 
t = time from start of flare. 


For constant t and increasing approach groundspeed (i.e., sink rate), the control law calls for a 
proportional increase in flare initiation height to achieve transient free flare initiation. 
Hence, in different wind conditions which give rise to different groundspeeds (i.e., different 
sink rates), the flare is initiated at different altitudes. Also, the time to touchdown from a 
given altitude is independent of speed. As a result, for varying groundspeeds, landing occurs 
at different positions along the runway. 


PRECEDING PAGE BLANK NOT FILMED 
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This problem has been investigated under the NASA TCV program (ref. 2) and one of the 
solutions was a “variable tau flare law.” In this flare law, the value for t (see equation 32) is 
varied as a function of groundspeed, i.e., 

V G0 


T — T. 


°V, 


(33) 


G 


where 


t = flare response time constant 
t 0 = nominal time constant 

V G o = nominal approach groundspeed 

V G = actual groundspeed 

Substituting Equation (33) into Equation (31) yields 

+ (A + A b ) = 0 

V G dt 

dh dh dx dh 

dt dx dt dx G 

Using Equation (35), then Equation (34) becomes 

T <Xao~r + (A + *b) = 0 
ax 

The solution to this equation is: x 

^oGo 

h(x) = (h 0 + h B )e - h B 


(34) 

(35) 


(36) 


(37) 


This solution shows that when Equation (36) is satisfied, then the resulting trajectory is inde- 
pendent of groundspeed, and only a function of distance. Also, flare may be initiated transient 
free at constant altitude, regardless of thrust speed. 

The implementation of a variable tau flare law was successfully test flown on the TCV air- 
plane and satisfied Category II landing requirements. Provided this concept could be imple- 
mented easily with the TECS, then it offered a proven solution in which the theory had been 
validated by flight test. 


6.2 VARIABLE TAU IMPLEMENTATION 

The equation used for implementation of the variable tau flare law can be derived by consider- 
ing Equation (31) from above: 


K h h + K „(h + h B ) = 0 

(38) 

r is defined as K h /K h , hence, 


T 'h + (b +^b) = 0 

(39) 


24 



Groundspeed variation can be compensated for by substituting in Equation (33), i.e., 


T = 


To 


YlLo 

v G 


( 40 ) 


Therefore, 


r 0 ^r°h + (h + h B ) = 0 (41) 

V G 

If the h signal is complemented with a height rate signal derived using radar altimeter mea- 
surements, then it is not necessary to utilize height directly. 

Hence, the control equation can be rewritten: 

V Go . 1 . h 1 . . 1 

To — h + - (h + h B ) = 0 or — + — — (h + h b ) - = 0 (42) 

V G 5 V G t 0 Vg 0 5 

However, since the flare time is inversely proportional to groundspeed, hb must be made 

proportional to groundspeed to assure equivalence of Equations (41) and (42). Implementation 
of this control law fits well in with the existing TEC architecture (see fig. 12). The justification 
for using h (derived from h R and hj) and not hj feedback is discussed in the following para- 
graph. 

To ensure a smooth, low-g maneuver, flare initiation should start at an altitude of approxi- 
mately 40 ft. A necessary feature of flare is to ensure a transient free switch over to flare from 
glideslope mode. For the nominal approach speed of 117 kn, height rate is approximately -10 
ft/s. To achieve transient free switch to flare, the gain ratio between the h-path and the inte- 
grated h-path (h-path), shown in Figure 12, should be 4. It should also be recognized that this 
gain ratio is the flare time constant (t 0 ) shown in Equation (42). 

The gain K FL = 1 /tqVqo in Figure 12, incorporates this gain factor r 0 . Compensation for 
groundspeed is achieved in the h-path as part of the basic TECS architecture. The bias term in 
the (h BIA g) path ensures the desired sink rate at touchdown. In addition to ensuring that the 
throttles retard during flare, a negative V CMD signal also is switched in at flare. The V CMD 
signal provides an additional elevator command that aids the rotation of the airplane during 
flare. 

The system as described above makes use of an h signal feedback. This derived signal is used 
in place of simply hj to enable the airplane to flare successfully on runways that are not 
horizontal. To cope with runways that slope, dip, or are humped, it is_necessary to utilize radar 
altimeter measurements complemented with hj signals to derive an h-signal in which the low- 
frequency components are referenced to the runway (fig. 12). 

6.3 RESULTS FOR VARIABLE TAU FLARE LAW 

The performance of the variable tau flare law was assessed using the flight simulator. It 
should be noted that h b was implemented as a constant. Figure 13 shows the flare law for 
nominal weight and normal approach velocity in zero wind. The control law has not been fully 
optimized to cope with all situations but for the present gain values shown in Figure 12, sink 
rate was 3.2 ft/s at touchdown, and touchdown displacement from runway threshold was 
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1,200 ft. In addition, the system was tested at minimum and maximum weight configurations 
using the appropriate approach speed for the weight (figs. 14 and 15). Sink rate at touchdown 
was 3 ft/s for maximum and minimum weight. 

These figures show that the flare dynamics are satisfactory for the various test cases, summa- 
rized in Thble 4. The flares are somewhat short for heavy weight and long for light weight. 
This will improve by use of the correct h b (proportional to groundspeed). The system has not 
been examined in headwinds or tailwinds with this configuration. 


7 able 4. Summary of Results for Variable Tau Flare Law 


RUN 

NO. 

GROSS WEIGHT 
(lb) 

V CAS 

(kn) 

HDOT AT 
TOUCHDOWN 

(ft/s) 

TOUCHDOWN 
ALONG RUNWAY 

SX(ft) 

TIME TO 
FLARE 

(s) 

1 

85,000 

117 

-3.2 

1,200 

6 


115,000 

139 

-3. 

1,150 

4.25 


70,000 

106 

-3. 

1,275 

7.25 


6.4 GAMMA-VDOT FLARE LAW 

Implementation of the variable tau flare law concept to the integrated autopilot/autothrottle 
system highlighted the basic requirements for flare: 

1) To rotate the FPA through an appropriate angle such that the sink rate is reduced to the 
touchdown requirement. This rotation should be independent of velocity (groundspeed). 

2) The throttles should be returned to idle during flare and rollout. 

The integrated autopilot/autothrottle system has been designed around a fixed inner loop 
configuration. All autopilot modes generate an FPA command (tcmd) an( f an acceleration 
command (V CMD ). It was considered that this inherent feature of the system could provide a 
means of obtaining flare. The inner loop configuration could remain intact and flare achieved 
by commanding a y signal to achieve the necessary rotation of the velocity vector of the 
airplane and obtain the desired sink rate at touchdown. In addition, the V command signal is 
used to retard the throttle and aid rotation of the airplane velocity vector. Tb arrive at a flare 
law in which the flare touchdown performance is independent of groundspeed, it was neces- 
sary to introduce this flare flightpath angle command through a lag. Since the FPA response 
must be quicker at higher groundspeeds to maintain the same flare trajectory, the lag time 
constant of tcmd must be scaled inversely proportional to groundspeed. The magnitude of the 
Tcmd signal is computed to give the required h at touchdown knowing V G . Also, the V CMD 
signal is scheduled proportional to V G to help minimize touchdown dispersion. 

As with the variable tau flare law, complementary filtering of hj and h R (radio altimeter 
measurements) was employed to provide for sloping runways. Figure 16 shows the implemen- 
tation of the gamma-Vdot flare law. 
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6.5 RESULTS FOR GAMMA-VDOT FLARE LAW 


To test the gamma-Vdot flare concept, an ACSL simulation was developed using the TCV 
B737 linearized aerodynamics for approach (nominal 120-kn EAS, flaps 40-deg). This simula- 
tion does not consider the effect of airplane geometry on touchdown distance (i.e., the vertical 
separation between center of gravity and the landing gear has been ignored). In addition, 
ground effects of the airplane have not been modeled. 

The gamma command time constant (t f ) was set to 0.75 sec and •h TD to -2.5 ft/s (fig. 16). This 
time constant (r F ) is combined with the basic -/-response time constant ( — 2.5 sec) to provide 
the desired flare response time constant ( — 4 sec) for the overall system. However, no attempt 
is made to force the system to accurately follow a specified trajectory. The results for no wind, 
15-, and 30-kn tailwinds are shown in Figures 17 through 20. Figure 17 shows height and 
height rate plotted against X, displacement along runway from the glideslope interception 
point. This point is 800 ft past runway threshold. Flare initiation occurred at an altitude of 40 
ft. Figure 18 is an enlargement of the touchdown area highlighting the dispersion in displace- 
ment (90 ft) and touchdown sink rate (0.2 ft/s). Figure 19 shows the lagged y response to 7 cmd 
and also shows the variation in groundspeed during flare. Figure 20 illustrates the throttle 
position during flare and also shows effect of flare on a and 8. These responses are all 
satisfactory. 

The sensitivity of the flare law to turbulence was investigated by adding longitudinal turbu- 
lence (5-ft/s rms). These results are summarized in Table 5. It can be seen that the standard 
deviation of the displacement for y-*V flare was 5% lower than the existing TCV airplane 
(standard deviation for H is 25% lower). 


Table 5. Effect of Horizontal Turbulence (U g - 5-ft/s rms) for y-V Flare Law 



y-V FLARE LAW 

EXISTING TCV AIRPLANE 


X G/S 

H 

X G/S 

H 


ft 

ft/s 

ft 

ft/s 

mean 

547.7 

-3.0 



s.d. 

199.9 

0.66 

211.2 

0.89 


• Results obtained from ensemble average of 15 runs (using ACSL simulation) 

• Corresponding TCV airplane results calculated from D6-37006; Lambregts and Hansen 
(ref. 2) 
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6.6 RESULTS OBTAINED USING THE FLIGHT SIMULATOR 


The y-V was also evaluated on the flight simulator. The results shown in Figures 21 and 22 
illustrate the performance for no wind and 32-kn headwind. Sink rate varied from ~ 3 ft/s for 
the no wind case to 2.4 ft/s for headwinds of 32 kn. The zero wind flare lasted 6.5 sec compared 
with 10 sec in a 32-kn headwind. A comparison of the y CMD an d 7 traces shows the effect of 
groundspeed on the flare time constant. 

6.7 GO-AROUND MODE 

The inner loop configuration of the TECS, which controls FPA and speed, allows for the imple- 
mentation of a simple but effective go-around mode. This implementation is shown in Figure 
23. The go-around is achieved by switching out the glideslope mode signal and switching in 
the desired positive FPA command. Simulation experience and results from pilot evaluation 
showed that a FPA of ~ 10 deg was satisfactory. This FPA provides a takeoff thrust go-around 
for weights greater than ~ 85,000 lb and reduced thrust go-around for lower airplane weights. 
For passenger comfort, normal acceleration levels may be reduced by lagging the FPA com- 
mand with a filter of time constant tq A . During go-around the system will continue to control 
speed to the selected command, without significant control coupling errors. 

The performance of go-around is shown in Figures 24 and 25 for two values of tq A . The 
airplane is on final approach down the glideslope and go-around is engaged at an altitude of 
100 ft. In Figure 24, tq A is small (0.2 sec), and g levels during go-around reach 0.5g. However, 
the height loss after engagement was only 24 ft. 

For large values of tq A (3.0 sec), maximum acceleration (AZCG) was lower at 0.3g, but height 
loss increased to 36 ft. The final value of tq a will be decided following further pilot evaluation 
of the system. 

In the simulation results, it can be seen that the system can just achieve the 10-deg FPA at the 
aerodynamic condition and weight below 1,000 ft. At about 1,000 ft, the throttle is at the 
forward limit and performance limitations are reached. Speed control has priority, hence the 
system maintains speed but the maximum attainable FPA (GAMMAD) drops below 10 deg. 
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7.0 VELOCITY VECTOR CONTROL COLUMN STEERING 


7.1 INTRODUCTION 

Velocity vector control column steering (vel-CCS) is a semiautomatic mode giving the pilot 
direct control over the inertial FPA angle yj through a column input. Column inputs are used 
in a rate command/hold system change the FPA reference command for the vel-CCS control 
laws and shape the system response to provide natural airplane control handling 
characteristics. 

The objective of the pitch vel-CCS mode is to reduce pilot workload when flying a desired 
flightpath. For example, in following an approach profile with defined FPAs but varying 
speeds, aircraft configuration, and winds, then a substantial reduction in pilot workload will 
result if the pilot can simply set up the desired inertial FPA path using column inputs and the 
automatic control system can track this command in spite of the disturbances. 

7.2 VELOCITY CCS CONTROL AND DISPLAY SYSTEM REQUIREMENTS 

Earlier work on developing vel-CCS for the NASA TCV airplane (ref. 3) led to the formulation 
of general control and display requirements for the vel-CCS: 

1) The response between the pilot’s input and the response of the display variable through 
which the pilot closes the inner control loop should be critically damped and exhibit a 
small enough lag to avoid pilot-induced oscillations. 

2) The y-response should be smooth, well dampened, and overshoot free. For steady state 
column deflection, the y-response should quickly establish the desired steady-state rate of 
change of FPA (fig. 26). 

3) Pitch attitude overshoot should not be excessive. Reference 3 established a lower limit of 
1 sec for the y-response lag to meet this requirement. 

7.3 DESIGN OF VELCCS SYSTEM 

The basic control concept of the TECS is to control FPA (y) and speed using control commands 
to the elevator and thrust. The vel-CCS develops a ycMD proportional to the integral of column 
deflection (fig. 27). The display uses both ycMD an d y. Use of ycMD to close the profile control 
loop bypasses the lag between command and response. The y CMD signal is used to drive the 
basic TECS algorithm without modification. 

One desirable feature was to maintain the same a z levels for the same column deflection as 
with the existing airplane. 

The TCV B737 airplane has a d e /g value of 15.89. This gives a column displacement (6 co j)/g of 
4.83 in at 120 kn, 40-deg flaps. 

Hence, 


b co( /'h = 0.15 in/ ft/s 1 


(43) 
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Since 


then 


7 = 


/?57.3 

V g 


deg/s 


(44) 


= — = 1.88 (45) 

5 co/ 0.15(202.2) in 

This ratio ylh col established the gain GAINV. In addition, a deadzone is incorporated in the 
output from the column. This was to ensure zero column signal to the control law when the 
column is in the neutral position. The results for this system are shown in Figures 28 and 29. 
Run 1 shows the basic response to the ramping gamma command. With the standard TECS 
gains, the 7 response lag is approximately 4 sec. Root locus analysis showed that the integra- 
tor gains in the thrust and elevator processing loops (KTI and KEI) could be doubled decreas- 
ing the response lag to approximately 2 sec. However, the resulting overshoot of the final 7 
value is undesirable. Furthermore, gain changes in the basic TECS algorithm for a specific 
mode are undesirable and will not reduce the 7 response lag to the desired level of 1 sec. 
Reducing the lag requires faster elevator response to rotate the airplane to the required FPA 
and tighter control of throttle to achieve a parallel 7 to 7 cmd response with a lag of ~ 1 sec. 

This was accomplished by feeding a proportional and derivative 7 cmd signal directly to the 
elevator and a proportional 7 cmd signal directly to the throttle command. 

Designing the gains for this configuration requires consideration of the physics of the prob- 
lem. For a response lag ( 7 ) of 1 sec and y c of 1 -deg/s, the steady-state FPA error (y e ) will be 
1 -deg (fig. 26). The command signals to the elevator and throttle can be calculated to achieve 
this steady-state error during column inputs, from which the system gains can be calculated. 

Considering the change in throttle position required to maintain 1 -deg/s rate of change of FPA 
with zero longitudinal acceleration, then 


AThrust = A7W 


where W = weight, 7 in radians/sec 


(46) 


and 


AThrust 

AT 


= 450 lb/ deg 


(AT = incremental throttle position) 


Hence, 


AT = 


A7W _ 1.85000 


= 3.3 deg/s 


450 (57. 3). 450 

For the existing system, the required y t of 1-deg contributes a throttle rate of only 

7 e #K TI deg/s = g 1.5 = 0.84 deg/s 


( 47 ) 


30 



Hence, a proportional path of 7 cmd t° the throttle should be added with a gain 


KTp2 = 5 7.3(3.3 0.84) = 4J? 

8 

During the steady-state y rate, the airplane altitude must change at the same rate as y. 
Therefore the rate of change of pitch command must equal y for the existing system. Pitch 
rate command $ c is 


* KEI 1 n „„ j / 

= 7 eg = 8 -7TT- = 0.662 deg/s 

a 57. 35 

Hence, a proportional path of y c must provide an additional signal of 
1-0.562 = 0.438 

and therefore, the gain of this proportional y c path (KGCP) should be: 

- 1 - KGCP g = 0.438 
57.3 

KGCP = —x 0.438 = 0.78 
8 

However, because of the gamma feedback signal of magnitude 1/57.3 rad/s, a gain component 
of unity is also required. Therefore, total gain required is 

KGCP = 1 + 0.78 = 1.78 
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In addition to the proportional y c path to elevator, a proportional y c signal path (KCPS) is used 
to provide a command signal to balance out pitch rate feedback during pitch maneuvers. 

From Figure 30: 


KCP g.Kfl.KEP = 1 K q 

57.3x 1.88 9 

Kq 57.3 x 1.88 4 x 57.3 x 1.88 

KCP = g K,KEP " 32.2 x 8.1 “ 1673 

However, these calculations were for steady-state conditions and adjustments were made to 
improve the transient response. The resulting gains were KTP2 = 3.9, KGCP = 2.75, KCP = 
2.0. For this system, a column deflection of 0.5 in. held for 8 sec gives a ramping Ycmd> which 
reaches a final y of 6.2-deg (fig. 31). The y response lag is reduced to 1 sec and the response is 
parallel to the 7 cmd- Overshoot of y is 3% and maximum error in CAS is 2 kn. Throttle 
response, velocity error, and attitude are shown in Figure 32. 

Figures 33 through 36 show the effect on performance of ± 50% variation in the significant 
gains of the vel-CCS system. As discussed above, the proportional and derivative gains (KGCP 
and KCP) affect the rate of rotation of the airplane to the 7cmd> the response lag, and to a 
certain extent overshoot (figs. 33 and 34). The proportional and integral gains to the throttle 
(KTI and KTP2) affect parallelism and final overshoot (figs. 35 and 36). 

The system was tested over a range of aerodynamic conditions from 120-kn EAS 1,500 ft, to 
310-kn EAS, 20,000 ft (figs. 37 through 40). The result showed that the gains associated with 
CCS control could remain fixed without degrading performance significantly. Overshoot in- 
creased to 4.8% in the high-speed case compared with 3% for the nominal case. Variation in 
parallelism was insignificant over the whole flight regime. 

7.4 DISPLAY CONCEPTS 

Requirements for the display have been met by the use of tcmd f° r use by the pilot to close the 
inner control loop. The actual FPA is also displayed for response monitoring. Display of both 
7CMD an( t 7 reduces the likelihood of pilot intervention during command tracking in 
turbulence. 

Previous work with CCS that examined the display requirements described the synthesis of a 
single display signal combining 7 cmd an d 7- This approach is feasible when the response lag 
of y is constant for all flight conditions and reduces display clutter. However, during periods 
with the column in detent, only actual y would be displayed and perturbations in turbulence 
may draw the pilot unnecessarily into the control loop. 
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8.0 PERFORMANCE ASSESSMENT 


The following paragraphs document the overall performance of the system as implemented on 
the Boeing Renton Flight Simulator [M-Cab]. 

8.1 ALTITUDE CONTROL 

For large maneuvers, an earlier altitude control mode was unsatisfactory in that it was impos- 
sible to maintain low g levels. In addition, h was limited restricting the maximum perfor- 
mance of the system. The problem with g levels is shown in Figure 41 with a +2,000-ft 
altitude change. It can be seen that maximum normal acceleration (AZ) reached -10 ft/s 2 and 
velocity error reached 5 kn. This problem was solved by the use of nonlinear elements de- 
signed to restrict acceleration levels and vary the start of exponential capture as a function of 
height rate. Figure 42 shows the effect of the nonlinear design in reducing normal accelera- 
tion levels. 

The performance of the system to a small height change (+300 ft) and a large change (+2,500 
ft) is shown in Figures 43 through 46 for low-speed and high-speed flight conditions. The small 
height change shows the linear performance of the system with 95% of the final value 
achieved by about 33 sec. Typically for small and large height changes, captures of the com- 
manded height are obtained without overshoot with maximum a 2 levels less than 0.1 g and 
velocity error less than 0.5 kn. At the extremes of the performance envelope (e.g., fig. 46), 
when commanding a lower height such that the throttles remain at the aft limit for a signifi- 
cant time, the nonlinearities of the engine spool up affect the throttle response giving rise to 
higher a z levels than desired (~0.15g). This peak a z lasts for only a very short time. This 
problem will be examined in the future. 

8.2 VELOCITY CONTROL 

The velocity hold system allows the airplane to control Mach, CAS, or groundspeed. The speed 
error (e.g., Mach or CAS) is converted to a TAS error which is processed through a gain and 
limiter to generate a V C md which drives the inner loops. The major system changes, different 
from that described in Reference 1, are the redesign of the gain and limiters so that for any 
combination of command inputs of speed and FPA, the system will satisfy the performance 
requirements discussed in Section 4. 

The performance results for the velocity control system are presented in Figures 47 through 
49. It can be seen that typically the height error does not exceed 10 ft. The linear performance 
(fig. 47) is constant over the performance envelope and the commanded velocity is captured 
with an exponential capture (r = 12.5 sec). 

The ability to control Mach is shown for a high-speed, high-altitude situation where control of 
Mach is mostly used (fig. 49). The control dynamics are identical to the CAS mode. 

The system is designed to automatically switch between CAS and Mach control during climb 
and vice versa during descent. This feature is discussed in Reference 1. 
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8.3 a CONTROL MODE 


The primary requirement of the a control mode is to prevent stall by holding the airplane to 
the safe a reference should the angle of attack reach this reference. This can occur when the 
pilot commands too low a velocity. The system described in Reference 1 worked successfully 
during approach at a fixed flap setting, but was not developed to cope with all flap settings 
and over the full flight regime. This has now been accomplished by varying the a reference as 
a function of flaps, varying the loop gain as a function of speed, and by linear filtering and the 
addition of nonlinearities to eliminate adverse throttle reaction. The mode was tested at differ- 
ent velocities, airplane weights, flap settings, and CG positions. Typical results are shown in 
Figures 50 and 51 for a limiting in response to commanding a low velocity. Figure 50 shows 
the response with maximum airplane weight, flaps fully extended, and aft CG position (31%). 
Figure 51 shows the result of minimum weight, 0-deg flaps and forward CG position (5%). 
Maximum height error recorded was 17 ft (fig. 51). 

8.4 FPA CONTROL MODE 

The performance of the FPA control mode is demonstrated at low speed (figs. 52 and 53). The 
linear performance for ±2.5-deg FPA command (fig. 52) shows overshoot free capture with 
approximately 2-sec time constant. Velocity error is less than 1 kn and vertical acceleration 
levels (a z ) were less than 0.1 g. 

For large FPA commands, a z levels are limited by the addition of a rate limit in the FPA error 
(y f ) path. This can clearly be seen in Figure 53 where the effect of the limit is shown in a z . 

In the case of large FPA changes, it is necessary to modify the control law configuration when 
throttle limiting occurs. This problem has been discussed in detail in Reference 1. However, 
briefly restating the problem: When the throttle limits, only the elevator remains available 
for control, and therefore only one variable can continue to be controlled and the preferred one 
is speed. The approach maintains the decoupled control concept whereby control command for 
one variable does not cause significant errors for the other controlled variable. This speed 
priority also avoids the danger of stalling and overspeeding. To this effect, the FPA error 
crossfeed to the elevator integrator is temporarily disconnected until the FPA command is 
reduced or the target altitude is captured. The response to large 7 cmd with throttle limiting is 
shown in Figure 53. In each case, the FPA settles at the maximum value commensurate with 
maintaining velocity. 


8.5 GLIDESLOPE CONTROL MODE 

The basic control law to capture the ILS glideslope has been discussed in Reference 1. This 
system has not been changed significantly. Performance is illustrated in Figures 54 and 55. In 
Figure 54, the engagement starts in altitude hold mode at 2,000-ft altitude approximately “2 
dots” under the glideslope. This corresponds to an initial glideslope error of 0.66 deg. The 
glideslope is captured in a overshoot free manner with neglible speed transient with the 
airplane following the glideslope to an altitude of 40 ft, at which point flare is initiated. In the 
noise free simulation run shown, the glideslope error is negligible. 

A flaps-up glideslope capture with incremental flap extension is presented in Figure 55. 
Again, the glideslope is captured overshoot free due to the fully adaptive FPA transition 
feature. This was achieved as a result of the concept of normalizing the glideslope error into 
the FPA command. This design approach also eliminates glideslope track errors due to decel- 
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eration. At the start of the run, a low velocity is commanded and the airplane flies at a ref 
(flaps 0-deg). Angle of attack (ALW) can be seen to approach a re f. As flaps are extended, 
velocity decreases and the angle of attack follows the new a ref at each new flap position. The 
glideslope error is negligible throughout this procedure. 

8.6 VERTICAL PATH MODE 

Vertical path mode allows the airplane to follow preprogrammed straight vertical paths in 
space. Switchover between vertical paths is arranged to occur in a transient free manner and 
the airplane captures each new vertical path in an exponential manner, with no overshoot. An 
example of this control mode is shown in Figure 56. The airplane has been preprogrammed to 
capture 5,000 ft, fly straight and level until it is over a waypoint, then climb at an FPA of 0.1 
rad to intercept a second waypoint at 7,500 ft. After a period of straight and level flight and a 
speed increase, the airplane passes over a third waypoint and descends at an FPA of 0.05 rad 
to reach the fourth waypoint at 5,000 ft. The flightpath and speed command executions are 
virtually perfect. 


8.7 CONSTANT ENERGY MANEUVERS 

The TECS has been designed on energy considerations. An airplane flying at constant height 
and velocity has an energy level consisting of the sum of the kinetic energy (KE) and potential 
energy (PE), where KE is associated with speed and PE with altitude. Certain maneuvers 
require a net energy change (e.g., an increase in velocity or height). This requires mainly 
control of thrust. In gusty conditions or when executing a double maneuver (i.e., decrease in 
height coupled with increase in velocity), the net energy change may be small but the system 
requires a redistribution of energy. The elevator is controlled to provide the means of redistrib- 
ution of KE and PE. 

This type of double maneuver in which the change in PE is equal to the change in KE is 
shown in Figure 57. The figure shows the command to decrease altitude by approximately 600 
ft is matched (in terms of energy) by the command to increase velocity 20 kn. As no net energy 
change is required, throttle motion is negligible and the elevator provides the primary means 
of redistributing energy. 

8.8 CLIMB AND ACCELERATE MANEUVERS 

A double maneuver requiring simultaneous climb and speed increase is illustrated in Figure 
58. The airplane climbs from 10,000 to 20,000 ft during which the speed changes from 250- to 
300-kn CAS. It can be seen that the airplane pulls maximum normal acceleration (a z ) of 0.1 g 
until the steady-state climb rate at maximum thrust is reached. Height rate peaks at 80 ft/s, 
then decreases slightly as altitude increases. On input of the commanded change in velocity, 
height rate decreases as a consequence of 50% of the available energy rate being used to 
increase velocity. Priority is given to velocity control and it is a simple matter to adjust the 
amount of energy rate used to execute the speed change (50% in this case). Once the required 
velocity is reached, then all available energy rate is converted to climb rate to reach the 
commanded altitude. The altitude is captured exponentially without overshoot. Throughout 
the maneuver, acceleration levels do not exceed 0.1 g. 
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8.9 DESCENT AND DECELERATION MANEUVERS 


The opposite maneuver to climb and acceleration is descent and deceleration. This is shown in 
Figure 59 for a 4,000-ft descent from 15,000 ft coupled with a 50-kn decrease in speed from 
300-kn CAS. In this case, the throttles retard to the aft limit, and the established sink rate 
depends on the drag configuration of the aircraft with the elevator continuing to control speed. 
When the command to decelerate is dialed in, velocity control has priority, height rate de- 
creases to zero, and 100% of available energy is directed to capturing speed. The use of 100% 
of available energy to capture speed was agreed on during pilot evaluation. However, this 
percentage can be varied easily with TECS to suit agreed flying procedures. In this situation, 
deceleration is constrained to prevent the airplane increasing altitude while velocity is being 
captured. Once the required velocity is reached, then sink rate increases again until the new 
altitude is captured. 


8.10 TURBULENCE AND WINDSHEAR 

The velocity control loop employs a complementary filter to provide a derived V feedback 
signal. The time constant of this filter largely dictates the final windshear and turbulence 
performance for the current system. The high-frequency component of V consists of the de- 
rived V signal obtained from 

AThrust . _. , AO , 

— — 7 g (see fig. 60). (48) 

Mass 

The first term represents longitudinal acceleration due to thrust and the second term is the 
component of gravity along the flightpath. 

The low frequency component of V consists of rate of change of airspeed. This derived V signal 
was developed to eliminate adverse pitch response to sudden gusts and undesired throttle 
responses due to flap and gear lowering. Variation of the filter break frequency (1/r) dictates 
the tradeoff between windshear and turbulence performance. A low frequency breakpoint 
severely attenuates the air mass desired V signal, thereby filtering turbulence and resulting 
in low throttle activity. However, as the high-frequency components have been attenuated, the 
response of the system to a windshear has been lagged, and as a consequence, the velocity 
error due to windshear is higher. Obviously, with a high-frequency breakpoint, the situation is 
reversed. 

Simulation runs were made with either windshear (1-kn/s headwind) or turbulence modeled 
using the Dryden wind spectrum. Moderate levels of turbulence were used (i.e., 5-ft/s rms 
longitudinal or horizontal turbulence, and 2. 5-ft/s rms vertical turbulence). The two selected 
aerodynamic conditions were: 1) low speed (125-kn CAS, 500 ft, flaps 40 deg, gear down), and 
2) high speed (325-kn CAS, 25,000 ft, flaps 0 deg, gear up). 

Figures 61 through 64 show standard deviation (SD) of velocity, height, throttle, and elevator 
against complementary filter time constant (ry) for longitudinal turbulence. Except for veloc- 
ity, the variables show a marked increase for increase in filter breakpoint (i.e., reduction in 
ry). Throttle shows high sensitivity at high speed, whereas the elevator shows increased activ- 
ity at low speed. In the case of the throttle within the engine control loop, the gain between 
throttle and thrust increases with increase in altitude. Therefore, the throttle shows increased 
sensitivity to variation in ry at high speed. The opposite situation exists for the elevator, 
where the gain between altitude command and elevator decreases with increase in speed. 
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In windshear (without turbulence) the maximum velocity error increased with increasing ry 
(lowering breakpoint), whereas height error decreased (figs. 65 and 66). Thus, improved speed 
control is possible at the expense of path tracking. 

Plotting the standard deviation of throttle position in turbulence against maximum velocity 
error due to windshear (fig. 67) shows clearly that a tradeoff can be made between the throttle 
activity and velocity error. At low speed, good windshear performance is essential to prevent 
stall and high throttle activity is the penalty. At the high-speed cruise condition, greater 
emphasis is placed on low throttle activity and windshear performance can be sacrificed. 

Ib achieve this aim, ry was varied as a function of altitude. For the high-speed case, it was 
decided to set r v at 20. For larger values of ry, throttle activity does not decrease significantly 
whereas velocity error due to windshear increases rapidly. 

For the low-speed situation, ry was set to 10 sec. Although it would have been desirable to 
have decreased ry still further to reduce velocity error in windshear, it was determined in 
simulations of landing approach that high elevator activity and rate limiting of the elevator 
made the performance unacceptable. Figure 68 illustrates this sensitivity of elevator to varia- 
tion of ry. 

Variation of complementary filter time constant ry has little effect on rms control activity in 
the presence of vertical turbulence. The effect of vertical turbulence was examined at low and 
high speed, the results of which are summarized in Table 6. 


Table 6. Effect of Vertical Turbulence (W g = 2.5-ft/s rms) 


LOW SPEED HIGH SPEED 

(HO=500 ft, CAS=125 kn) (H0=25,000 ft, Mach=0.74) 


Filter breakpoint (ry) Filter breakpoint (ry) 

5 10 20 5 10 20 


Elevator (deg) 

0.33 

0.31 

0.33 

0.08 

0.07 

0.07 

Throttle 
Position (deg) 

1.61 

1.55 

1.68 

1.14 

1.11 

1.05 

Vel error (kn) 

0.48 

0.59 

0.88 

0.10 

0.14 

0.14 

H error (ft) 

4.02 

3.63 

3.96 

3.14 

2.88 

2.67 
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9.0 CONCLUSIONS 


The integrated longitudinal autopilot/autothrottle based on the total energy control system 
(TECS) concept has been developed and evaluated in detail. The feasibility of the TECS design 
has been demonstrated by meeting all design objectives with a system that is significantly 
simpler than existing conventional designs. In addition, the system design features numerous 
system advantages and performance improvements over conventional designs: 

1) Use of a single generalized flightpath and speed control algorithm for all modes and 
flight conditions allows elimination of most autothrottle and flight management com- 
puter control laws for the longitudinal axis. 

2) The number of computer interfaces is reduced substantially compared to conventional 
designs. 

3) Except for the elevator and throttle innerloop control, the design is completely general 
and transferable to other airplanes without further development or adjustment. 

4) The system operation and performance is consistant for all modes and flight conditions, 
reducing implementation checkout and certification flight testing requirements. 

5) Speed and flightpath control are decoupled for all modes and flight conditions, such that 
speed errors due to flightpath maneuvers are less than 1 kn and path deviations due to 
speed maneuvers are less than 20 ft. 

6) Energy exchange maneuvers are executed through the elevator, without significant 
throttle response. 

7) Vertical maneuvers are g-limited, except for the go-around mode. 

8) Flightpath intercept maneuvers are fully adaptive to the flight condition, resulting in 
overshoot free captures at all times. 

9) No significant flightpath deviations are incurred due to longitudinal accelerations/ 
decelerations. 

10) No significant speed deviations are incurred during climb or descent. 

11) During maneuvers in which the thrust limits, elevator priority is given to speed control 
over flightpath control, preventing overspeed or underspeed. 

12) A flare mode using the total energy control concept was developed. This mode commands 
an incremental FPA proportional to altitude change. A simultaneous deceleration com- 
mand is used to aid the pitch rotation and retard the throttles. 

13) A vel-CCS mode was incorporated in the TEC system concept, using gains tailored for fly- 
by-wire implementation. Performance equivalent to that achieved on the NASA TCV 
B-737 aircraft was shown. 

14) Performance in windshear and turbulence has been optimized. 

PRECEDING PAGE BLANK NOT FILMED 
39 (Page 40 is blank) 


10.0 REFERENCES 


1. A. A. Lambregts, “Vertical Flight Path and Speed Control Autopilot Design Using Total 
Energy Principles,” AIAA Paper 83-2239CP, August 1983. 

2. A. A. Lambregts and J. F. Creedon, “Development and Flight Evaluation of Automatic 
Flare Laws With Improved Touchdown Dispersion,” AIAA Guidance and Control Confer- 
ence, Paper No. AIAA-80-1757, August 1980. 

3. “ACSL Advanced Continuous Simulation Language User Guide and Reference Manual,” 
Mitchel and Gautier Associates, 2nd edition, 1975. 


PRECEDING PAGE BLANK NOT FILMED 


41 





Throttle EPRIDL 



EPRIDL 

Figure 2. Simplified EPR Control Loop 


EPR 



38 deg 


Figure 3. Simplified JT8D-9 Engine Model 
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Figure 11. Spool Up and Spool Down Characteristics 
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Figure 12. Variable t Implementation 
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Figure 13. Variable rFlare Law: 85,000 lb, V C/ 


54 


VC AS (kn) 



120 

110 

100 


-THCM (deg) 




THETAD (deg) 



= 117 kn, CG = 0.20 








0.00 1.50 3.00 4.50 6.00 7.50 9.00 10.50 12.00 13.50 15.00 

TIME (SEC) 


LEFT AXIS RIGHT AXIS 


Figure 14. Variable rFlare Law: 115,000 lb, V CAS = 119 kn, CG = 0.2 
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Figure 15. Variable rFlare Law: 70,000 lb, V CAS =106 kn, CG = 0.20 
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Figure 16. Implementation of Gamma-VDOT Flare Law 
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Figure 20. Theta, Throttle Demand, and Alpha Against Time 
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Figure 22. Flare Results: 32-kn Headwind 
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Figure 23. Gamma-VDOT Flare Law With Go-Around Mode 
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Figure 30. Velocity Vector-Control Column Steering System 
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Figure 33. Variation of Gain KGCP 
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Figure 34. Variation of Gain KCP 
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Figure 35. Variation of KTI 
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Figure 36. Variation of KTP2 


76 




□ Run 1 

<D Run 2 


EAS Alt 
120 kn, 1500 ft 
150 kn, 5000 ft 


THRDEM THETA versus T 




79 





I 


□ Run 3 
O Run 4 


EAS Alt 
200 kn, 10,000 ft 
310 kn, 20,000 ft 


THETA (deg) THRDEM THETA versus T 










3300 

3150 

3000 

2050 

2700 







Figure 43. Altitude Control Mode: 120 kn, 3,000 ft, Ah = +300 
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Figure 44. Altitude Control Mode: 120 kn, 3,000 ft, Ah = +2,500 
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Figure 45. Altitude Control Mode: 325 kn, 25,000 ft, Ah = +300 ft 
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Figure 46. Altitude Control Mode: 325 kn, 25,000 ft, Ah = -2,500 ft 
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Figure 47. CAS Control Mode: 140 kn, 500 ft, FL = 40 deg, AV = +10 kn 
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Figure 48. CAS Control Mode: 250 kn, 25,000 ft, FL = 0 deg, AV = +50 kn 
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Figure 49. Mach Control Mode: M 0.65, 25,000 ft, AM - +0.15 
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Figure 50. a Control Mode: 150 kn, 5,000 ft, 105,000 lb (31% CG, FL = 40 deg) 
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Figure 51. a Control Mode: 200 kn, 5,000 ft, 70,000 lb, 5% CG, FL = 0 deg 
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Figure 52. FPA Control Mode: 120 kn, 3,000 ft, FL = 40 deg, Ay = +2.5 deg 
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Figure 53. FPA Control Mode: 120 kn, 3,000 ft, Ay = +15 deg 
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Glideslope error 
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Figure 54. Glideslope Control Mode: H 0 = 2,000 ft, FL = 40 deg 
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Figure 55. Glideslope Control Mode With Flap Extension and Touchdown 
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Figure 56. Vertical Path Mode 
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Figure 57. Double Maneuver (Ae = 0) 
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Figure 58. Climb and Accelerate Maneuver 
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Figure 59. Descent and Deceleration Maneuvers 
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Figure 62. SD Height Against 7y 
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Figure 63. SD Throttle Position Against 
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Figure 64. SD Elevator Activity Against ry 
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Maximum velocity error due to 1-kn/s windshear (kn) 
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Figure 65. Maximum Velocity Error Due to Windshear Against t v 
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Figure 66. Maximum Height Error Due to 1-kn/s Windshear Against Ty 
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Figure 67. Maximum Velocity Error Due to 1-kn/s Windshear (kn) 
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